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1. Introduction 


Pattern formation in systems of large size has attracted much research interest in recent 
years. Especially the fact that many aspects, at least of one-dimensional systems or 
of quasi one-dimensional patterns, can be described by reduced equations of motion 
has allowed for linking quite different fields of physics (cf. JI] and references therein). 
To some extent the approach strongly parallels the normal form calculations in low 
dimensional dynamical systems 0, The reduced equations for simple instabilities, 
i. e. the Ginzburg-Landau equation is well established and its derivation can even 
be found in textbooks ||. However, at least from the general point of view, less is 
known if additional constraints are imposed on the instability, that means if higher 
order codimension bifurcations are considered. 

We are here concerned with instabilities in spatially one dimensional systems 
where a single mode becomes unstable with respect to a wavenumber q c due to an 
eigenvalue zero in the spectrum. Such a situation occurs generically in inversion 
symmetric situations and we henceforth consider such systems. Let A (k,/i) denote the 
corresponding critical eigenvalue in dependence on the wavenumber k and the system 
parameters //. It obeys 


A(g c ,/r) = 0, d k X(k,fi) t =0, d 2 k \ 


k=q c 


k=q c 


< 0 


( 1 ) 


These conditions trace a codimension one set, i. e. a hypersurface in the parameter 
space on which the instability occurs. It is well established, even from a rigorous point 
of view (D, that the dynamics near such an instability is governed by a slowly varying 
envelope which obeys a real Ginzburg-Landau equation. Whether the instability is sub- 
or supercritical, i. e. whether the amplitude saturates in the vicinity of the instability, 
depends on the sign of the cubic term. The transition from sub- to supercriticality, 
i. e. the change of the sign, leads to a codimension two bifurcationf. It is of course 
contained in the bifurcation set determined by eqs. ([[]). Such instabilities, which for 
example are relevant in the hydrodynamic context (cf. [5] for a recent reference) are at 
the centre of interest of our contribution. 

From pure symmetry considerations the structure of the reduced amplitude equation 
may be fixed, taking the translation and inversion symmetry into account 


d T A = rjA +DdlA + f\A\ 2 A + s\A\ 4 A + ig\A\ 2 dzA + idA 2 dtA* . (2) 

However, these considerations do not tell us whether such an equation is valid at all, 
and how the coefficients depend on the actual parameters of the underlying equations 
of motion. We present the complete derivation of the amplitude equation (j2j) starting 


f Sometimes such a bifurcation point is called a tricritical point, in order to distinguish it from a 
codimension two bifurcation caused by the degeneracy of two distinct modes. 
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from a general equation of motion, even if the method is in principle well established in 
the hydrodynamic context. But our approach is purely algebraic and has the advantage, 
that the results can be applied immediately to quite different physical situations. As a 
by-product we remark that for the similar hard-mode case a comparable approach fails, 
in contrast to statements in the literature. Finally we will dwell on some properties of 
eq.(|2|), since a complete discussion is difficult to find in the literature, despite the fact 
that related results from different points of view can be found quite frequently SBBB- 

2. Derivation of the amplitude equation 

2.1. Notation 

We suppose that the basic equation of motion for the iV-component real held t ) is 
cast into the form 

d t $ = C ! & + W[$;/r] (3) 

such that the trivial translation invariant stationary state is given by $ = 0 f. For the 
linear operator, which determines the instability of this state, we allow for an expression 
as general as possible, i. e. 

= . (4) 

a 

Using plane wave solutions uexp(ikx) the eigenvalue problem is completely determined 
by the N x N matrices 

= ( 5 ) 

a 

according to 

= A(fc,/£)«*(//) . (6) 

We denote by A the eigenvalue branch with maximal real part, which obeys eqs.(||) at 
the instability. In our case the inversion symmetry guarantees that all quantities are 
real valued. 

For the nonlinear part we employ an expansion in powers of the held according to 

A/(£; /i) = A/" 2 (£; fi) + TV's OF; /r) + A/^l; /i) + A/" 5 (£; //) + • ■ ■ (7) 

where 

A40E;h) = (8) 

a,f3 

ms.; H) = £ S, d?s, m) (9) 

a,13 ,7 


f In particular we concentrate on situations where boundary conditions play no significant role, so that 
we can consider formally systems of infinite extent. 
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denote the most general expressions of second and third order, with vector valued bi- 
and trilincarforms C_^ a ^ and . Written in components they read for example 

( Q[ { “ p) {u, n}) m = 51 c ( ^ kl u k vi . (10) 

k,l 

The contributions of order four and five are understood in the same way using the 
notation and foA/fofo for the corresponding multilincarformsf. In the subsequent 

analysis it is necessary to evaluate the nonlinearities if plane waves are inserted for the 
field. To be specific only waves with the multiples of the critical wavenumber q c will 
occur. In such a case all the nonlinearities are expressed in terms of the abbreviations 
(cf. eq.(|)) 

C mn (u,v]fj,) :=Y,( [m< lc) a (mqc)^C^ ) {u,v} (11) 

afS 

D lmn (u,v,w-,i±) :^^{Uq c ) a (imq c )^(mq c yD^ 7) {u,v,w} (12) 

<787 

E klmn (u,v,w,x;[j,) := ^ {ikq c ) a {ilq c ) 0 (imq c y(inq c ) s E^ S) {u,v,w,x} (13) 

—jkimn(M-i V-i W-i2Li yj f) ■= J2 (y9c)“(ifc?c)^(ii?c) 7 (i^?c) i (iwgc) e Zj* /J7fc) {M,lZ,m, J/}(!4) 

aP'ySe 

which are frequently used in what follows. 

2.2. Weakly nonlinear analysis 

Suppose that at /i = // a degenerated soft-mode instability occurs. Let w c exp(k/ c :r) 
denote the marginally stable mode, i. e. u c = Uq c (F c ) is the nulleigenvector of the matrix 
(||) at n — fi and k = q c . In the vicinity of this parameter value, i. e. for 

j£ = + £ f/ 1 ^ + £ 2 /r (2 - ) + • • • (15) 

the solution of the full equation is expanded as 

<t(x, t) = £ 1/2 $i + e %2 + £ 3/2 $3 + £ 2 $4 + £ 5/2 $ 5 H- (16) 

$i = u c exp(iq c x)A(T 1 ,T 2 , • • •, 6, 6, ■ ■ ■) + c - c ( 17 ) 

where £ denotes a dimensionless smallness parameter. As usual the dynamics of the 
complex valued amplitude A, which possesses a slowly varying space-time dependence 
on the scales r k = e k t, = e k x, will be determined by the secular conditions of the 
expansion. However, before we proceed let us comment on the choice of the expansion. 
The actual expansion parameter is given by e 1 ^ 2 and for completeness the slow scales 
e x l 2 t and e 1 ^ 2 x should also be taken into account. But as will become obvious from the 

f In addition the symmetry properties C^ a ^{u, y} = u}, D^ a ^{u,v,w} = 

qjhccS) 

[w, u, v, } = ..., ... are employed in what follows. 
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following considerations these scales drop and do not contribute. The same conclusion 
holds for the expression (|15|), where // 1,a ) act as unfolding parameters. In addition, 
since /d 1 ) will be confined along the codimension one set of soft-mode bifurcations, 
the real unfolding of the bifurcation occurs at the order 0(£ 2 ). Hence the scaling 
of the amplitude in eq.(16) coincides with the scaling in the corresponding spatially 
homogeneous situation. Nevertheless, we stress again that the inclusion of all terms of 
orders 0(e k ^ 2 ) yields the same results as presented below. 

If one inserts the expansion (0) into the definitions (|) and (j7) one obtains 

C tL = C + £jC {1) + e 2 £ (2) + ... (18) 

and 


= JV*(£) + eA/?’® + ... 


(19) 


Here the terms of higher order contain the parameters // 1,2 h In order to simplify 
the notation we do not introduce a different symbol for the contributions of order 
zero, but just skip the argument fi. Even this convention introduces a slight abuse of 
notation, it is henceforth understood that the corresponding expressions are evaluated 
at n = n , e. g. C := . Analogous expansions hold for quantities like eq.(||) or 

f 11) f 14). In particular the matrix L(k) = L(k: u ) determines the critical mode and 

i«(A) = (£ (1 >9 e )i(fc;rtl H = iv 


The following steps are now like for the usual codimension one case and can in 
principle be found in textbooks. If one inserts the expansion (16) into the equation 
of motion (|H), takes eqs.([L8|) and ([L9|) into account and performs the derivatives with 
respect to all scales, then one obtains at each order in e k A an inhomogeneous linear 
equation determining 


exp (inq c x)w n . (20) 

n 

The inhomogeneous part typically contains Fourier modes with integer multiples of the 
critical mode, and the slow scales are just considered as fixed parameters. If v c denotes 
the left nullcigenvector of the critical mode, i. e. v*L.(q c ) = 0, then the condition that 
the solution of eq.(^Oj) does not become secular reads 


fecla-i) = o 


( 21 ) 


Here the brackets denote the usual scalar product. Now the solution, discarding 
exponentially decaying transients, reads 

c h, = - J2 L(nq c )~ 1 w n exp(mq c x) 

n^± 1 

~ ( K L~ 1 w n=1 exp([q c x)+u £ Z(r ll ...^ 1 ,...)exp(iq c x) + c.c.^ (22) 
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where the coefficient Z of the solution of the homogeneous equation may depend on the 
slower scales. LA 1 denotes the inverse on the subspace omitting the critical mode, i. e. 

£%(*)= UQc)^ 1 = 1-IMcXeJ, 47'Mc = 0, v* c L- 1 = 0 . (23) 

The reason that we reiterate this scheme here is twofold. On the one hand we would like 
to give the reader the chance, within the amount of formalism to relax with a passage, 
which he or she of course knows quite well. On the other hand the textbooks mentioned 
above usually stop with such general considerations or specialise to certain conditions 
which may not be shared by the model under consideration. Here we will continue with 
the most general equation of motion, even if we have to proceed to the fifth order. We 
demonstrate that the explicit evaluation is not at all horrible within a suitable notation. 

At the order 0(e 1//2 ) eq. (|20[) just yields the eigenvalue equation for the critical mode 
(cf. eq.(^)). At order O(^) the quadratic nonlinearity contributes nonresonant Fourier 
modes ±2g c , 0 to the inhomogeneous part (cf. eq.([A3|)), so that no nontrivial secular 
condition occurs. The solution (22) reads in this case 

$ 2 = 2T 20 |A|- + (r 22 exp(i2 q c x)A 2 + u c exp(iq c x)B + c.c.) (24) 

where the amplitude B of the homogeneous solution depends on the slower scales and 
the abbreviations ( |A4| ) are introduced. 

At the third order 0(£ 3 ^ 2 ) the quadratic and cubic nonlinearity contribute as well 
as the linear operators £ and d t} if the derivatives with respect to the slow scales are 
performed. The complete inhomogeneous part of eq.(22) is for convenience given in 
eq.(A6). The secular condition (^T|) yields 

o = -(Vcl Uc)d ri A + (v_c\L (l \q c ) Uc)A + p 3 A\A\ 2 (25) 

if we equate the coefficient of <%A to zero with the help of eq.([l]). In the usual 
codimension one case we would now end up with the Ginzburg-Landau equation. Here 
however we require that the coefficient of the cubic term vanishesf 

:= (WISGsid^,^) +4Ci 0 (M c ,r 20 ) + 3Z2m(u c ,u c ,u c )) = 0 . (26) 

Together with the condition O this equation determines the codimension two 
bifurcation manifold. Since the secular condition (^5|) has to yield a finite and 
nonvanishing solution, we finally have to require that both of the two remaining terms 
vanish separately. Hence we are left with 

d n A = 0 (27) 

and 

/x (1) <9 i£ A(^ c , /x) U= Hc = 0 (28) 

f Indices with an overbar denote negative values n := —n. 


if we take into account, that the matrix element in eq.(^5|) can be expressed in terms of 
a directional derivative owing to the definitions (|l^) and (p~ 8 [) . The condition (|28|) has 
a simple geometrical interpretation (cf. fig .[[]). If we take the total derivative of eqs.([I|) 



Figure 1 . Diagrammatic view of the parameter space in the vicinity of the degenerated 
soft-mode bifurcation point /i c - The solid line indicates the soft-mode bifurcation 
manifold and /F 1 ’ 2 ) the unfolding vectors (cf. eq.(|l5|)). 


with respect to n along a direction in the codimension one bifurcation manifold, i. e. we 
take the dependence of the critical wavenumber on /i into account, we are exactly left 
with eq.(^). Hence the secular condition fixes the parameter variation /b 1 ) a t order 
O(^) in such a way that only variations within the soft-mode instability manifold are 
permitted. The full parameter unfolding is obtained at higher order. By the way we 
remark that similar considerations exclude the spatial scale e x ^x from the perturbation 
expansion. For the solution we now get the result 


$3 = 


L 3 iA\A \-- i7 31 <%^ + 7 b ^A exp(i q c x) + T 33 exp(i3g c a;)H 3 


-L31 


2 T 20 AB- 


+ 2T 22 exp(i2 q c x)AB + w c exp(i q c x)C + c.c. 


(29) 


using the abbreviations (A 8 )-(A 10 ). 

At the order 0(k 2 ) also the parameter dependence of the nonlinearities (cf. eq.(|19j)) 
contributes. The inhomogeneous part, which is given in eq.(Al2) up to Fourier modes 
2 q c yields for the secular condition (| 21 |), taking the secular condition (|28|) of the preceding 
order into account 


0 = — {Vc\Uc)d T1 B + p 3 2\A\ 2 B + A 2 B* 


By virtue of the higher order codimension condition ( [261) we are left with 


d T1 B = 0 


(30) 

(31) 


For the solution at this order we have, if we restrict to Fourier modes ±2 q c , 0 which will 
become resonant at the next order 







+ exp(i2g c x) r 42 A 2 | A|~ + 2v^ 2 Ad^A + 2r 22 AC + t^ 2 A 2 + H 22 -S 2 + c.c.'j + .. .(32.) 

The coefficients are given by eqs.(Al3)-(Al8). 

We now ping in all results to compute the secular condition at the order 0(£ 5 ^ 2 ) 
and obtain, taking eqs. ([!]), (|27|), 


(33) 


), and (pTf ) into account 

(Uc\Uc) —d T2 A + r]A + Dd^A + i cd^A + rA|A | 2 + sA|A | 4 + ig\A\ 2 d^A + i dA 2 d^A 
+p 3 \2\A\ 2 C + A 2 C* + A*B 2 + 2A|5| 2 1 - (v c \uf)d T1 C = 0 . 


Thanks to the higher order codimension condition (|26|) the nonlinear terms which conple 
the different amplitudes vanish. Furthermore the last summand, being solely dependent 
on the scale r 4 has to vanish too, in order to avoid a secular contribution. If we introduce 
A := exp[ic£i/( 2 D)\A to eliminate the linear derivative term, we are left with the closed 
amplitude equation (0), where 


g := q + c 2 /(4D), r := r + c(g — d)/(2D) 


(34) 


It is worth to mention that our formalised approach has enabled us to incorporate the 
higher order codimension condition at all steps in the perturbation expansion. 


2.3. Coefficients 


In addition we have obtained the general microscopic expressions for the coefficients. 
We use the notation introduced in section 0 and the abbreviations of the appendix. 
The linear unfolding parameter reads 


(Mc\Uc)v ■= (vc\k a) + {Vc\L i2) (q c )Uc) 

= (vcllLc) [l/ 2 (h ( 1 ) ^) 2 A(g c ,/i)|^ c + (p ( 2 ) dff)X(q c ,g) 


(35) 


where the last expression follows from the definitions flT5|) and ([18]) straightforwardly. 
Hence the linear unfolding contains a contribution from the curvature of the bifurcation 
manifold and one from the transversal intersection. 

The diffusion constant is given by 


{VcK)D := —l/2(v c \ff\q c )u c ) - (vfffiqc)^) = -l/2(v c \u c )dl\(k, gj\ k=qc . (36) 

For the linear derivative term we have obtained 


(VcK)c ■= {He I -■*^ (1) (?c)7^ 1 - C L { 1 ) y{qc)Uc ) 

= ~ (Vc\Ue)d k (g {1) d^)\(k : g)\ k=qc ^ = ^ . (37) 

In view of the relations (Jj) the derivative can be expressed also in terms of the change 
of the critical wavenumber along the soft-mode bifurcation manifold. 
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The cubic unfolding coefficient reads 

(vcK) r ■= {vc\L (1) (qJLsi + 2(Z 2 i(r 22 ,7^) + 2C 2 I (75 2 ,m c ) 

+ +4C 10 (u c ,tJ 0 ) +GD 11 i{^ v y^,y^) + 3D u i(y^,u ci ^ 1 ) 

+ 2C$ (T 22 , u c ) + 4£Zio (m c ,E 20 ) + 3^!(m c , M c , M c )) . (38) 

If we use here the representations (A10), (A20), and (A21) as well as the higher order 
codimension relation (|26|) the expression simplifies to 

{Vc\Uc)r = (h (i) ^)P3(g c ;h)U c =M ( 39 ) 

if we define the object on the right hand side by eq.(^ 6 j) but evaluated with the full 
parameter dependent nonlinearities (11), (12) and eigenvectors (cf. eq.(||)). 

The evaluation of the coefficient of the quintic term yields 

{lLc\Uc) s := (lfc|2hZio(M c ,£40) + 2C 2 i(r 22 , T 3 i) + 2C 2 j(T 42 ,m c ) + 2(23 2 (r 33 , T 22 ) 

+ ^iC^oQZsi, T 20 ) + 6 ^nj(T 31 ,Mo u^) + 3J2ni(u c , u c , T 31 ) 

+ 3 -D 3 jj(££33, ILci ILc) + 6i2 2 i 2 (T 22 , M,,, T 22 ) + 12D 100 (u c ,T 20 ,r 20 ) 

+ I2£ 20l (r 22 ,r 20 ,u c ) + 4^ 1112 (« c , McjHc? r 22 ) + i2A 21 n(r 22 ,u c , « c , m c ) 

+ 24i? 110l (u c ,u c ,r 20 ,u c ) + lOZmnfe,^,^,^,^)) • (40) 

This expression is a genuine term of the fifth order and cannot be reduced further. 

For the normal derivative term the coefficient reads 

(Hcl Uc)g ■= {Uc\ - 2^ / (g c ) (9c) T 31 - 4C( 10 (u c , 7 ^ o ) + 4C 2I (7“ 2 ,m c ) 

— 4^ 10 (t^ 1 , T 20 ) — 4C 01 (T 20 ,u c ) — 4C( 10 (u c , T 20 ) — 4C 2 i(T 22 , Me) 

~ Me,Me) - OZZmtMcMcMc)) (41) 

whereas for the odd derivative term we obtain 


(lLc\lLc)d ( 22 cI — A / (q , c)E 3 1 + 4 Ci 0 (u c , 21 (—22; 7 ^) 

— 4C 0 i(r 20 , Me) — 2 Ci 2 (u c , T 22 ) + 3 ^ni(u c ,u c , 7 ^ 1 ) — SDin^c, u c , u c )) . (42) 

All coefficients are real valued, since the constituents are real owing to the symmetry 
of the underlying system. The coefficients of the linear terms can be expressed in terms 
of the spectrum. In addition, the cubic unfolding is obtained as a formal parameter 
derivative of the cubic coefficient of the ordinary Ginzburg Landau equation. One 
should note that this contribution and the linear derivative term are both caused 
by the parameter variation along the codimension one bifurcation manifold and are 
easily missed if the parameters are not unfolded according to eq. (pA|). The mentioned 
properties are not passed to the amplitude equation (pjj, since the coefficients are 
renormalised by eqs.([l4|). For the remaining coefficients no simple interpretation seems 
to be available. 
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3. Properties of the amplitude equation 

Partial discussions of eq.(j2j) from different points of view can be found in the literature 
Here we focus on those results which have in our opinion consequences for the 
behaviour near the codimension two bifurcation point. 

First of all the coefficients D and s have to be positive respectively negative in order 
to yield a bounded solution. We confine the subsequent analysis to this case. Hence 
these coefficients can be incorporated in the length scale as well as the magnitude of 
A and an additional parameter can be eliminated by a rescaling of the time. However, 
since no real simplification is achieved we discuss the unsealed equation directly. 


3.1. Potential case 


In the absence of the odd derivative term, d — 0, eq.(0) admits a potential L = f £dx 
decreasing in time, with density 


t := - g\A\ 2 + D\d^A\ 2 - r/2\A\ 4 - s/3|H | 6 - ig/4\A\ 2 (a* 8^A - Ad^A*) 
= - fj\A \ 2 + D d^A + ig/(4D)\A\ 2 A 2 - f/2\A\ 4 - [s/3 + g 2 /(16D)} \A\ 


(43) 


The potential is de fini te for s < —3 g 2 / (16.D), so that in some sense every solution tends 
to a time independent state. However, if the inequality is violated, e. g. if g is to large, 
the solutions may diverge. The potential property seems to be destroyed, if an odd 
derivative term is present. 


3.2. Bifurcation scenario 

The trivial state A = 0 of the amplitude equation is stable if fj < 0. Beyond this 
threshold time independent plane waves emerge from this solution. For the existence of 
these solutions A = a K exp(ffi£) we obtain from eq. (Q) the condition 

0 = fj — Dk 2 + (f — AiA)\a K \ 2 T s|a K | 4 . (44) 

The quantity A := g — d completely incorporates the dependence on the normal and 
the odd derivative term. 

Consider for the moment an arbitrary but fixed wavenumber k. Eq. ( |44| ) determines 
the bifurcations of the corresponding plane wave. It is evident that at 

fj = Dk 2 (45) 

a wave emerges from the trivial solution, and it is generated in parameter space on that 
side of the bifurcation set where the inequality 

(fj — Dn 2 )(r — An) < 0 


(46) 
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is valid. Hence this peculiar bifurcation changes from sub- to supercritical behaviour 
at f = A k, i. e. at 

fj = Dr 2 / A 2 . (47) 


Now we are considering larger amplitudes a.k and concentrate on the case, where the 
waves exhibit a saddle node bifurcation. If we rewrite eq.(fH|) in the form 

r ~ A ^ 2 tj-Dk 2 ' ^ 

we immediately recognise that a saddle node bifurcation occurs at 

9 / . (f — A k,\ 2 

V — Dk — (s) ^ _ 2g J , (49) 

provided that the inequality 

(f- Ak)/(-2s) > 0 (50) 



holds. The content of eqs.(f45p, (f46[ ), (j47|), ([49P , and (|5Q| ) is summarised in fig.|] for a 
particular wavenumber k. The region of existence of a plane wave solution, which is 
bounded by eq. (f49|) extends beyond the stability region of the trivial state. 

We now have to perform the analysis presented above for every wavenumber k. The 
region where plane wave solutions exist is given by the union of the regions described 
above. Hence its boundary is determined by the envelope of the curves (fh^) for all 
wavenumbers. It is easily computed as (cf. fig.Q) 


rj = -Dr 2 /[4D{-s)-A 2 }, (f> 0), if |A| < A c := 2 sjD(-s) .(51) 


If |A| approaches A c the parabola ( 0 ) degenerates with the negative 77 -axis. For 
| A| > A c no boundary exists at all, so that for every 77 < 0 there exist plane wave 
solutions. 

I 11 summary, the scenario for |A| < A c resembles the sub-supercritical transition in 
low dimensional dynamical systems, where a saddle-node bifurcation line is typically 
born. But in the extended case this behaviour is destroyed at |A| = A c , which can be 
viewed as a codimension three bifurcation point. 

Yet we have not claimed the stability of the plane wave solutions. A linear stability 
analysis according to A — exp(i/c£)(l + 8(5) yields the linear equation 

d T 8(3 = fj8/3 + D(—k 2 8(5 + i2 nd^8f3 + 8^5 (3) + f\a K \ 2 (8(5* + 2 8/3) + s|a K | 4 (25/3* + 3 8(3) 

— (g — d)n\a K \ 2 (28(3 + 8(3*) + ig\ot K \ 2 8^8(3 + id\a K \ 2 8^8 (3* . (52) 


Splitting into real and imaginary parts and and taking the fixed point equation 
into account the corresponding real two dimensional system reads 


8 -r 



2 (r — AK,)\a K \ 2 + 4s|a K | 4 + Ddj (—2 Dk, — A| a 


(2-Dtc + (g + d)\a K \ 2 )d^ 


D81 



.(53) 
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Figure 2. Sketch of a partial bifurcation diagram of the amplitude equation (|2j) in the 
r—fj plane. The thick broken line indicates the instability of the trivial state. The thin 
broken line represents the bifurcation of plane waves from the trivial state for a wave 
number k with A k > 0 (cf. eq.p5|)). The pitchfork like inserts indicate whether the 
bifurcations are sub- or supercritical and the thin dotted line gives the transition for 
all wavenumbers (cf. eq.(Hd)). The thin solid line marks the saddle node bifurcation 
for the plane wave (cf. eq.(|49|)). The corresponding features for the wavenumber — k 
are displayed in grey. Finally the thick solid line marks the envelope of saddle node 
bifurcation lines for all wavenumbers (cf. eq.(|5l|)). 


Analysing the stability in terms of plane waves exp(ku£) yields a two-dimensional 
eigenvalue problem, where the trace and the determinant of the corresponding matrix 
read 


2 — Du 2 + (f — Ak)\ a K | 2 + 2s \a K 4 

(54) 

u? D 2 u 2 — 2D (f — Ak) \a K 2 — 4Ds\a K 4 
(2Dhi + (C 7 + d)\a K \ 2 j ( 2 Dhi + A|o ; re | 2 j 

(55) 


Stability requires Tr^ < 0 and Det^ > 0 for all wavenumbers u. Owing to the simple 
dependence on u the condition on the trace results in 

\a K \ 2 > (f — Ak)/(-2s) , (56) 

which is of course valid if the right hand side is negative. Whenever the right hand side 
is positive and the parameters are such that plane wave solutions are possible, i. e. we 
are beyond the saddle node bifurcation line (cf. eqs.(ff5|), (joD|)), then eq.(^) tells us, that 
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the solution with the larger amplitude obeys the constraint (|56|), whereas the solution 
with the smaller amplitude is unstable. Hence we are left with checking the condition 
on the determinant which results in 


4D(-s) \a K 


\ot K 


(r — Ak)/(- 2s) — [2Dk + (g + d)\a K \ 2 ^j (2 Dk + A|a K | 2 ) >0 -(57) 


Due to this condition it is evident, that the stability properties depend on the parameters 
g and d separately. A complete discussion of the stability using eq. (|57]) is straightforward 
but tedious. We do not intend to discuss the full implications of this inequality, but just 
concentrate on a neighbourhood of the envelope (0) in order to study whether stable 
solutions are generated at this bifurcation line. For that purpose we fix the wavenumber 
to hi = — A f / (A 2 — A 2 ) which is just the value for which the saddle node bifurcation line 
touches the envelope (cf. £g.|2| and eqs.(f7g|), ([TU|), and m and expand the left hand 
side of eq.(|57l) for fj values slightly beyond the envelope (cf. eq.|49|). We obtain to the 
leading order the result 


[4. D(-s) - 2gA\ \a h 




— (r — Ak)/(—2 s) 


> 0 


(58) 


For the solution with the larger amplitude (cf. eq. (]56|) ) the condition is satisfied provided 
that gA < 2 D(—s) holds. Then a stable plane wave occurs at the envelope. 


4. Conclusion 

We have presented the complete and systematic derivation of the amplitude equation, 
which governs the transition from sub- to supercritical soft-mode instabilities. Within 
our approach the general expression for the coefficients and especially their dependence 
on the system parameters has been obtained. Although these formulas look a little bit 
lengthy, one should keep in mind that they can be applied to almost every physical 
situation, and that their evaluation is straightforward in concrete cases. 

From the principal point of view it is worth to mention, that the amplitude equation 
(@) can be derived consistently at all. Such a feature is far from being obvious. To 
emphasise this point consider the corresponding hard-mode case, i. e. an instability 
at q c = 0 with a nonvanishing frequency. A superficial inspection would suggest that 
the whole derivation goes along the same lines with minor modifications. But if we 
follow the approach of section [|, we are left at the third order with eq.(^). Since 
now all expressions are complex valued but the higher order codimension condition 
requires a vanishing real part only, the secular condition becomes a nonlinear equation. 
Of course it can be easily integrated to yield the time dependence on the scale T\ as 
A = Aexp[iIm( / 03 )|A| 2 ri], Here the constant of integration A depends on the slower 
scales. If one uses this representation in the subsequent orders, then the derivatives with 
respect to spatial coordinates yield linearly in T\ increasing terms, since the exponent 
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is space dependent. This feature invalidates the systematic derivation although a semi- 
quantitative approach has been proposed (cf. the discussion in [|TQ| ). 

The origin of such difficulties lies in high frequency components which contribute 
to the secular conditions in low orders and cause an uncontrolled mixing of different 
time scales. Similar phenomena are well known in the problem of counterpropagating 
waves, where a formal derivation is still possible and a nonlocal coupling in the 
amplitude equation is generated |TT| . For the degenerated hard-mode instability 
we expect finally similar effects but further investigations are needed. Nevertheless 
these considerations emphasise again the necessity of careful derivations of amplitude 
equations to supplement phenomenological approaches. 

Concerning the behaviour beyond the sub-supercritical soft-mode instability we 
stress that without an odd derivative term a potential system occurs. Hence that 
kind of term is responsible for a persistent time evolution beyond the threshold. In 
addition, we mention that the difference A between the normal and the odd derivative 
term determines the domain of existence of spatially periodic patterns in the vicinity of 
the threshold. These properties again show that the behaviour beyond the instability 
depends crucially on the actual numerical values of the coefficients. 


Appendix A. Inhomogeneous part 


For convenience we list in this appendix the inhomogeneous parts, which occur in each 
stage of the derivation of the amplitude equation. For notational simplicity the same 
label is assigned to the components at each order and we use an overbar to denote 
negative indices n = —n. All the abbreviations which we introduce are real valued. 

If we insert eq.(16) into the equation of motion (|3j) and observe the expansion (p~8|) 
we obtain at the order 0(^ 1 ^ 2 ) 


0 = A(<?c)w c exp(zg c a;)A + c.c. 


(Al) 


since the spatial derivatives act on the plane wave only. This condition is fulfilled by 
virtue of the eigenvalue equation for the critical mode ([]). 

At the order O(^) one obtains a contribution from the temporal and spatial 
derivatives acting on $ 2 and one contribution from the quadratic nonlinearity (8) at 
p = /i with the derivatives acting only on the plane waves. Taking the abbreviation 
(11) into account the result reads 


9i«2 = 


a0 


= C %2 + £ZnCMc;Hc)|A| 2 + £Zn(HcHc) exp(2 iq c x)A 2 + c.c. 
Hence in the notation of eq.(^0|) the nonvanishing Fourier components read 
Mo = 2Z/(0)r.2oI> M 2 = ~ A(2q , c )r 22 A , 


(A2) 

(A3) 
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where the abbreviations 


£.20 ~A( 0 ) 1 Q.ll(lLcilLc)i 


£ 2 -A(2g c ) 1 £11 £0 Me) 


(A4) 


have been used, and the obvious relation w n = uf n should be noted. The solution of the 
linear equation (A 2 ), discarding exponentially decaying transients, is given by eq.(24). 

At the order 0(e 3 ^ 2 ) the time derivative and the linear operator act plainly on $ 3 . 
In addition, these derivatives give a contribution when acting on the slow scales of dq. 
For the nonlinearities now the quadratic and the cubic terms at /i = /i contribute 


d t § 3 = £$ 3 + £ (1) £ + (££) [1] - ££ 

+ E 2 {82*3, af + £ DW d^} 


(A5) 


a/3 


afi 7 


Here the notation (£$ 3 )^ means that the derivatives have to be evaluated at the first 
order in e. All other derivatives with respect to x are understood at fixed values of 
Taking the relation (<9“ exp(iq , c x)A)[ 1 l = a(iq c ) a ^ 1 exp(iq c x)d^A into account the third 
contribution on the right hand side is expressed in terms of the derivative of the matrix 
(H) with respect to k. If we evaluate the nonlinear contribution with the help of the 
solution ( |24| ) of the preceding order and recast all contributions into the form (^0j) we 
get 


Ml 0 — — 2^(0)r 20 AH 

ULi = ~ Ucd T1 A + ^(q^Uc-A — iL!(qc)Uc9^A 

+ [2C 2 i(r 2 2, Me) + 4C 10 (m c , r 20 ) + 3Z) n j(M c , Mc> Me)] |A| 2 A 

w _2 = — 2^{2q c )T 22 AB 

m = ~ A(3g c )T 33 A 3 (A6) 


with the abbreviation 


£3 := -L 1 (3q c )[2C 21 (T 22 ,u c ) + D lll ( 2 ^,v^, 2 i c )] (A7) 

Here L/(k) denotes the derivative with respect to k. The nonsecular solution discarding 
transients is given by eq.(29) where the abbreviations 

£1 := — L. c 1 [2 £ 2 i(£ 2 , ^c) + 4 (ZioCm c ,£o) + 3J2 n j (u C ’,UciUc)\ (A8) 

TSi := = (1 - \yJ)(Vc\)dkUk(n)\k=q c (A9) 

7^ := -L^L^iq^Uc = (1 - |Mc)(£cI)(^ (1) ^m)^ c (m)U=/£ c ■ ( Al °) 

have been introduced. 

Proceeding to the order 0(£ 2 ) one has to observe that in addition the parameter 
dependence of the quadratic term has to be taken into account (cf. eq. fll5D ). Denoting 
the corresponding directional derivative by | M=(U the equation 
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reads 

a* = csu + £ (1, #2 + (£4 2 ) [11 - a„i 2 

+ 5 : [c<“«R“i 2 , affe} + 2 c<-»{^3, a^} 

a0 

+ ^ 3 D<“W{aj 42 .^*i,a 2 fc}+E -(An) 

a /37 ol(3^5 

Inserting the solutions of the preceding orders and performing the derivatives we obtain 
for the Fourier modes up to wavenumber 2q c 

m = L{ 0) [-2r 20 (AC* + A"C + |B| 2 ) - 2r4 (A%/1* - A*%A) 

- L0I2II 4 - 2EJ2II 2 ] 

w_i — ~ Uc9 T1 B + ^■ 1 \q c )u c B — [^(q^Ucd^B 


+ [20.21 (£. 22 ,He) + 4Ci 0 (m c , r 20 ) + 3D.ni(Uc,24c,Me)] 2\A\' 2 B + A~B* 
m2 — - L( 2 q c ) 7^2 a 2 + 2 r 22 AC + 2 i 7 “ 2 A + r^kFA 3 + r 22 ir 

with the abbreviations 

£40 := — L.(Q) 1 [ 2 Cii(r 31 , Uc) + 2 C_ 1 i(u c , r 31 ) + 2C 22 (r 22 , r 22 ) + 4 C 00 (r 20 , r 20 ) 
+ 3 Z 2 112 (m c , Mo E 22 ) + 3 £> 2 n(r 22 , Uc, Uc) + 12i2ioi(M c , £20 j lie) 

+ 6Ei 1 ii(y c ,u c ,u c ,u c )] 

T^o := -^(O )- 1 [^( 0 )r 2O + C 11 (u c ,£ 1 )-C ll (u c ,u c ) 

A ■= -i(O) -1 [GiiC&.sU + fZii(s*,;&) +i (1) (0) Eao+di’fecSc) 

£ 42 — £(2?c) 1 | 2 C n (r 3 ii u0 + 2 C 3 i(r 3 3 l u,.) + 4C 20 (r 22 .r 20 ) 

+ 6Z2 110 (Mc) Uci r 20 ) + 6 i 2 2 ll(r 22 , Me, 2ic) + 6£liiu(Uc, Ulci Ulci Me)] 

7^2 := ~M 2 9c) _1 - £ii(w c ,M c ) 

lL := -i( 2 ? c )-‘ [2e 11 (a c ,75 1 )+£ (1) (29 c )r 22 +cSi ) (tt ; ,3fc)" 

and 


(A12) 


(A13) 

(A14) 

(A15) 

(A16) 

(A17) 

(A18) 

(A19) 


C mn (u,v) :=iJ2a(imq c )“ 1 (inj c )' 3 C ( “' 5| {a,ji} . 

a /3 

The nonsecular solution of eq.(All) up to Fourier modes ±2g c is then given by eq.(32). 
We remark that, if the dehnitions (A4) are understood in terms of the full parameter 
dependent quantities and eigenvectors (cf. eqs.(|]), (11), and (|6|)), then the abbreviations 
(A15) and (A18) obey, taking relation (A10) into account 


A= 

A = (k 1, 7)E2 2 


- r 20 (<iij(a (1, 7)ii7A)li!-i> + c.c. 

- 2E 22 (ial(ii (1) ap)!i IC (/i)U-».) • 


(A20) 

(A21) 
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Finally at the order 0(£ 5 ^ 2 ) we end up with 

5 = £$5 + £ (1) $3 + £ (2) $i + (^ 3 ) [1] - 5 T1 $3 + (£$i) [2] - <9 t 2 $1 + (£( 1) $ 1 )W 

+ £ [2C' (a/3) {«9“^4,5f$i} + 2C (a/3) {<¥$ 3 , df$ 2 } 

a/3 

+ 2C' ( “ /3) {(«9“^ 2 ) [11 , + 2C' (a/3) {(«9“ti) [1] , df$ 2 } + 2C' (q/3)(1) {^$ 2 , 

+ y [3D ( “' 3i »{s."43,sf4.,a:ii} + 3D<“«{s»2,^t2,;?*i} 

a/37 

a/3^8 

+ E . (A22) 

a/3^f8e 

As before dAAX 1 ) denotes the directional derivative and (£$i)^ indicates that spatial 
derivatives have to be evaluated at order 0(e 2 ). Inserting the previous orders and 
performing the derivatives the Fourier mode uy of the inhomogeneous part of eq.(A22) 
is evaluated and yields after some algebra the secular condition (33), taking the 
abbreviation 

j2imn(u,v,Ul) '■= i a(i/g c )" _1 (h«g c ) /3 (ing c ) 7 iZ (a/37) {^^lh} (A23) 

a/37 

into account. 
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